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Abstract 

We propose a partitioning of the set of unlabelled, connected cubic graphs into two 
disjoint subsets named genes and descendants, where the cardinality of the descendants 
is much larger than that of the genes. The key distinction between the two subsets is 
the presence of special edge cut sets, called crackers, in the descendants. We show that 
every descendant can be created by starting from a finite set of genes, and introducing 
the required crackers by special breeding operations. We prove that it is always possible 
to identify genes that can be used to generate any given descendant, and provide inverse 
operations that enable their reconstruction. A number of interesting properties of genes 
may be inherited by the descendant, and we therefore propose a natural algorithm that 
decomposes a descendant into its ancestor genes. We conjecture that each descendant can 
only be generated by starting with a unique set of ancestor genes. The latter is supported 
by numerical experiments. 



1 Introduction 



In this paper, we study the family S := S(N) of undirected, unlabelled, connected cubic graphs 
on N — 2m vertices, where m is a positive integer greater or equal to 2. We recall that a cubic 
graph has exactly three edges incident on every vertex. In our study we will be concerned with 
the problem of generating complex cubic graphs by appropriate compositions of simpler cubic 
graphs in a manner somewhat analogous to genetic "breeding" . This will be achieved with the 
help of six "breeding operations" . Cubic graphs that cannot be seen as resulting from such 
operations will be called "genes", and all other cubic graphs will be called "descendants". Of 
course, descendants constitute a majority of cubic graphs. We will prove that every descendant 
graph can be created from a finite set of genes. Since a lot of the structure in a descendant is 
inherited from the genes, many properties of those genes are inherited as well. This motivates 
the introduction of "inverse operations" that, ultimately, decompose a given descendant into a 
set of "ancestor genes". Such a decomposition could subsequently be used with other graph 
theory algorithms to improve solving time. 

Undirected cubic graphs have been extensively studied in the literature (e.g., see Harary 69 , 
|Balaban and Harary 71| [Holton and Shechan 93 ). In particular, there are now programs that 
enumerate all instances (up to an automorphism) of cubic graphs on N vertices (e.g., see 
|Meringer 99| ) . A significant line of research concerns generation of cubic graphs with the help 
of an exponential generating function (e.g., see [Brinkmann 96] . |Meringer 99| , |Royle 87] ). To 
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the extent that, in this paper, we generate descendant cubic graphs, the present contribution is 
conceptually, but not methodologically, related to the preceding. 



The construction of ancestors and descendants of cubic graphs could be seen as being related to 
the work of McKay and Royle McK ay and Royle 86| . However, a feature of our approach is that 
all ancestors and descendants remain in the class of cubic graphs, and collectively constitute the 
entire set of connected cubic graphs. This enables us to study those properties of cubic graphs 
that may be "inherited" from simpler cubic graphs. Furthermore, the construction methods 
are quite different from those in |McKay and Royle 86| . The work of Batagelj Batagelj 136] is 
also similar to ours in that it is concerned with the generation of complex cubic graphs from 
simpler cubic graphs, where the former maintain the properties of the latter. In fact, some of 
the generating rules given in Batagelj 86 are analogous to some of the breeding operations in 
this paper. However, the two approaches differ in that Batagelj 's approach involves starting with 
a single cubic graph and replacing particular structures within the graph sequentially. In our 
approach, the majority of development occurs by combining multiple cubic graphs together, with 
the resulting graph being more complex than both parents. 



Importantly, preliminary numerical experimentation revealed a surprising "unique ancestors 
property" that is either very common or, as we believe, universal. Namely, we conjecture that: 
Given any descendant cubic graph Gd there exists a unique set of genes from which Gd can be 
derived by a finite sequence of the six breeding operations. Of course, there may well be a mul- 
titude of sequences of breeding operations that lead from the ancestor genes to Gd- However, 
if the above conjecture holds, then all possible decomposition pathways of Gd into its ancestor 
genes will ultimately lead to a unique set of ancestor genes. 



2 Preliminaries 



At the conclusion of this manuscript, we conjecture that any given descendant can be obtained 
via our approach through a particular set of ancestor genes, and provide numerical evidence 
supporting this conjecture. If true, this conjecture, combined with the proof of existence of 
ancestor genes for any descendant, implies that our approach permits a unique decomposition 
for every cubic graph. 

A standard concept in graph theory is that of edge connectivity. In the simple case in which 
removing a single, specific edge would disconnect the graph, that edge is called a bridge. Extend- 
ing this notion a set of edges EC constitutes an edge cut if their removal disconnects the graph. 
The smallest cardinality of an edge cut in a given graph G is defined to be the edge- connectivity 
of that graph. Of course, cubic graphs can be k— edge-connected for only three values of k: 1, 2 
or 3. 

We note, however, that in some 3— edge-connected cubic graphs (e.g., the famous Petersen graph, 
see Figure [2]) the removal of any minimal edge cut set isolates a single vertex. Arguably, such a 
partition of vertices is in a sense degenerate, and prevents a more refined classification of cubic 
graphs. To address this problem and achieve a finer classification we introduce a special class of 
edge cut sets that we name crackers. The latter are defined as follows. 

Definition 2.1. An edge cut set of a cubic graph T consisting of k edges is a k- cracker if no 
two edges are adjacent in the sense of being incident on the same vertex, and no proper subset 
of these edges disconnects the graph. 
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Lemma 2.2. The removal of a cracker from a connected graph results in exactly two disjoint 
components. 

Proof. Since any cracker is an edge cut set, the removal of a cracker must result in at least two 
disjoint components. Assume there exists a cracker C — {ei, . . . , e^} whose removal results in 
more than two disjoint components. We can think of the removal of C as a series of individual 
edge removals, for each edge in C. It is clear that the removal of a single edge cannot result 
in more than one additional disjoint component. Therefore, after removing edges e\, . . . , e^—i, 
but before removing edge e^, there must already be at least two disjoint components. However, 
this implies that a proper subset of C also disconnects the graph, which violates Definition 12.11 
Therefore, C is not a cracker, and the initial assumption is false. □ 

Of course, any cracker is not only an edge cut set, but in fact a cyclic edge cut set. This is clear 
because its removal disconnects the graph into two connected subgraphs, each of which must 
contain at least three vertices. Since each remaining vertex has degree at least two (as only 
non-adjacent edges were removed), these connected subgraphs must contain cycles. It is then 
clear that a minimal cyclic edge cut set in a given cubic graph is a cracker of minimal size for 
that graph. However, larger cyclic edge cut sets may contain adjacent edges, and therefore not 
all cyclic edge cut sets are crackers. 

We recall that the girth, g, of a graph is the length of the shortest (nontrivial) cycle in the 
graph. In many cases, a nontrivial cycle determining a cubic graph's girth automatically defines 
a g-cracker made up of edges that are not in the cycle but which have one vertex on the cycle. 
However, many graphs have crackers of size less than g. For instance, it is easy to see that in 
the graph given in Figure [TJ g = 3 but the edges e\ and ei form a 2-cracker and no 1-cracker 
exists in this graph. 




Figure 1: An 8- vertex cubic graph with girth 3 and the two cracker {ei, e 2 } 

We note also that there are only two connected cubic graphs, on 4 and 6 vertices respectively, 
that contain no crackers at all (see Figure |3] and the discussion in Section 

Remark. Note that, since the minimal cyclic edge cut set in a cubic graph is a cracker, it is clear 
that for any given cubic graph, the cyclic edge connectivity is equal to the size of the smallest 
cracker in that graph. For the sake of simplifying the notation, we will refer to a cyclically 
fc-connected graph as a C/c-connected graph. The class of all Cfc-connected cubic graphs on N 
vertices will be denoted by Sk(N), or simply by Sk, when the number of vertices is fixed. 

We note that the famous Petersen graph is C5-connected in the above sense (see Figure [5]) , as 
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the edges connecting the "inner-star" to the outer boundary form one of a number of 5— crackers, 
and no smaller crackers exist in this graph. 




Figure 2: Petersen graph which is C5-connected 



3 Motivation 

It follows immediately from Remark [2] that (for fixed N > 8) the class C of connected cubic 
graphs can be partitioned as 

M 

s=\Js l , 

i=l 

N m . . N 

where M < — . The choice of upper bound is conservative because at most — non-adjacent 

N 

edges can be chosen in any graph of size N, but in reality it is likely that far fewer than — 
partitions will be required for any given N. For instance, when N — 20, M — 6. 

Definition 3.1. If a cubic graph F is C/c-connected for k > 4, we call it a gene. Otherwise, we 
call r a descendant. 

The reasoning behind the choice of names gene and descendant is made clear in Section 21 where 
we demonstrate that any descendant can be obtained from a set of genes, through the use of 
prescribed breeding operations that introduce crackers into a descendant. Experiments have 
shown that genes are far less numerous than descendants. 

It was mentioned earlier that two cubic graphs, namely the 4- vertex gene T| and the 6- vertex 
gene Tg, contain no crackers at all. These two graphs can be seen in Figure [3] The following 
lemma proves that every other cubic graph contains at least one cracker, and that the size of the 
smallest cracker is bounded above by the girth of the graph. 
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Figure 3: Cubic graphs T| and fg, which contain no crackers 



Lemma 3.2. Except for and Tg, all connected cubic graphs contain at least one cracker of 
size no more than the girth g of the graph. 

Proof. The cases of T| and Tg can be confirmed by inspection. There is one other cubic graph 
containing 6 vertices, with girth 3, which contains a 3-cracker, as displayed in Figure 2J So the 
lemma is true for N < 8. 

For any cubic graph containing 8 or more vertices, it was proved in Lou et al [Lou et al. 01] that 
there exists at least one cyclic edge cut set of size g. The smallest cyclic edge cut set in a cubic 
graph is a cracker, so it is clear that the smallest cracker can be of size no bigger than g. □ 




Figure 4: A 6-vertex cubic graph, with girth 3, which contains a 3-cracker comprising the high- 
lighted edges 

In Section[5l it is proved that any descendant graph can be obtained from a set of genes. It is our 
hope that descendants inherit many of their properties from the genes used to construct them. If 
so, any analysis of such a descendant could be reduced to the problem of analysing the component 
genes, which are often much smaller than the descendant. The subsequent investigation of graph 
theoretic properties in sets of genes is a natural topic for future research. 

One such graph theoretic property of interest is that of Hamiltonicity, that is, the property of 
containing a simple cycle of length equal to the number of vertices in the graph. We observe 
that non-Hamiltonian genes are extremely rare. Even excluding the (trivially non-Hamiltonian) 
bridge graphs, non-Hamiltonian descendants constitute a large majority of the remaining non- 
Hamiltonian graphs; see Table [TJ The second column of that table, labelled by NHi, lists the 
percentages of bridge graphs relative to the total cardinality of non-Hamiltonian graphtQ denoted 
by NH. The third column labelled by NH2+ := NH \ NHi, lists the percentages of graphs two 
or more cyclically edge connected relative to the cardinality of NH. The fourth column labelled 
by NH4+, lists the percentages of graphs four or more cyclically edge connected relative to the 

1 For a more complete study of the prevalence of cubic bridge graphs relative to the total set of cubic non- 
Hamiltonian graphs, see Filar et al |Filar et al.~10] . 
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cardinality of NH. Finally, the fifth column labelled by NH4+/NH2+, lists the percentages of 
graphs four or more cyclically edge connected relative to the cardinality of all non-bridge, non- 
Hamiltonian, graphs in NH 2 + . We shall define graphs in NH44. as mutants, a name that properly 
reflects their exceptionality. For instance, we note from the fourth column of Table [T] that with 
N = 18 only 0.12 of one percent are mutants, which corresponds to two (out of 1666) non- 
Hamiltonian cubic graphs on 18 vertices. These two mutants are the famous Blanusa Snarks 
[Blanusa 46j . 

In fact, it is no coincidence that the Blanusa Snarks appear as mutants in this framework. In 
Read and Wilson plead and Wilson 98j , the definition of an irreducible Snark is given as a cubic 
graph with edge chromatic number of 4, girth 5 or more, and not containing three edges whose 
deletion results in a disconnected graph, each of whose components is nontrivial. This final 
condition, along with the well known fact that all Snarks are non-Hamiltonian, is akin to our 
definition of a mutant. Therefore, the set of all mutants is a superset of the set of all irreducible 
Snarks. However, some non-Snark mutants do exist, and therefore have either girth 4 or an 
edge chromatic number of 3 (or both). In particular, the BH-Mutant, displayed in Figure [S] is 
the smallest non-Snark mutant, and has an edge chromatic number of 3. There are 16 further 
non-Snark mutants of size 22, one of which is the Zircon-Mutant, also displayed in FigureJS] The 
Zircon-Mutant also has an edge chromatic number of 3. 





NHi 


NH 2+ 


NH 4+ 


NH 4+ /NH 2+ 


10 Vertex 


50.00% 


50.00% 


50.00% 


100% 


12 Vertex 


80.00% 


20.00% 


0% 


0% 


14 Vertex 


82.86% 


17.14% 


0% 


0% 


16 Vertex 


84.93% 


15.07% 


0% 


0% 


18 Vertex 


86.13% 


13.86% 


0.12% 


0.86% 


20 Vertex 


87.40% 


12.60% 


0.05% 


0.38% 


22 Vertex 


88.59% 


11.41% 


0.02% 


0.21% 



Table 1: Distribution of Non-Hamiltonian (NH) Cfc-connected graphs 



4 Breeding and parthenogenic operations 

From Definition 13. 1[ it is clear that genes do not contain any 1-crackers, 2-crackers, or 3-crackers. 
Collectively, we refer to 1-crackers, 2-crackers and 3-crackers as cubic crackers. As a corollary, 
descendants must contain at least one cubic cracker. It then seems plausible that we might be 
able to construct any given descendant by combining two or more cubic graphs together in such 
a fashion as to create the cubic crackers present in that descendant graph. Since there are three 
different types of cubic crackers, we define three breeding operations that map two cubic graphs 
to a single descendant by inserting a cubic cracker between them in such a fashion as to retain 
cubicity. In such a case, we say that the descendant has been obtained by breeding. We refer to 
the original two cubic graphs as the parents of the descendant graph, and likewise the descendant 
graph is the child of the two parents. 

Note that the following operations are defined only for cubic graphs. Although they work for 
disconnected cubic graphs, in this manuscript we are interested only in connected cubic graphs, 
and make the assumption that all input graphs are indeed connected and cubic. 







Figure 5: BH-Mutant which is a girth 5, C4-connected non-Hamiltonian non-Snark mutant 
containing 20 vertices, and Zircon-Mutant which is a girth 5, C5-connected non-Hamiltonian 
non-Snark mutant containing 22 vertices. The colours identify Tait Colorings 



4.1 Breeding operations 

Definition 4.1. A type 1 breeding operation is a function B± defined on the tuple (ri, r 2 , e±, e-z), 
where Y\ — (Vi,Ei) and T 2 = (V 2 ,E 2 ) are cubic graphs, and furthermore, ei = (a, b) G E\ and 
e 2 = (c,d) G E 2 . This function maps such a tuple onto another tuple (Tr),e) as follows 



where = (Vd, Ed) and {e = (vi, v 2 )} G -Ed- The new set of vertices is Vn = V\ U V2 U «i U 
v 2 . The new set of edges is i?D = (-Ei\ei) U {E 2 \e 2 ) U {(a, ui), (6, Ui), (c, u 2 ), (d, U2), ^2)}- 

Note that a type 1 breeding operation always outputs a bridge graph (that is, a Cl-connected 
graph) . See Figure [5] for an illustration. 

Definition 4.2. A type 2 breeding operation is a function B 2 defined on the tuple {T\, T 2 , e\, e 2 ), 
where Y\ = (V%,Ei) and T 2 — (V 2 ,E 2 ) are cubic graphs, and furthermore, e\ = (a,b) G E\ and 
62 = (c, e?) G £2 and neither edge is a 1-cracker. This function maps such a tuple onto another 
tuple (r.D,e3,e4) as follows 



Z3i(ri,r 2 ,ei,e 2 ) 



(rz),e), 



£ 2 (ri,r 2 ,ei,e2) 



(r_o, e 3 , e4), 



7 




Figure 6: Graphs Ti, r 2 and as described in Definition 14.11 



where Tjj = (Vd,Ed) and {e^ = (a, c),64 = (b, d)} G Ejj. The new set of vertices is Vr> — 
Vi U V 2 . The new set of edges is E D = (E{\e{) U (-B 2 \e 2 ) U {(a, c), (6, d)}. 




Figure 7: Graphs Ti, r 2 and as described in Definition 14.21 



Clearly Yd contains the 2-cracker {(a, c), (6, d)}. See Figure[7]for an illustration. Note also that 
a type 2 breeding operation always creates a 2-edge-connected descendant, unless either of Ti or 
r 2 is 1-edge-connected (in which case, Yd is also 1-edge-connected). 

Definition 4.3. A type 3 breeding operation is a function B3 defined on the tuple (ri, r 2 , V\, t> 2 , a, b, 
c,d,e,f), where Ti = (Vi,E\) and T 2 = (V^-E^) are cubic graphs, and furthermore, v\ G V\ is 
incident to vertices a , b and c G V\ and w 2 G V 2 is incident to vertices d, e and / G V%- None 
of the edges adjacent to V\ or w 2 are 1-crackers. This function maps such a tuple onto another 
tuple (r^j, ei, e 2 , ea) as follows 

S 3 (ri,r 2 ,wi,w 2 ,a, 6, c,d, e, /) = (r D , ei, e 2 , e 3 ), 

where = (Vd,Ed), and also, {ei = (a, rf),e 2 = (6, e),e3 = (c, /)} G Ed)- The new set of 
vertices is Vn = (Vi\ui) U (V2\w 2 )- The new set of edges is Ed = (Ei\{(vi,a), («i,6), («i,c)}) U 
(£ 2 \{K d), (U2, e), («2, /)}) U {(a, d), (6, e), (c, /)}. 

See Figure [8] for an illustration of type 3 breeding. 



4.2 Parthenogenic operations 

In addition to the preceding three breeding operations, we also define three parthenogenic oper- 
ations. These are operations that map a single descendant to a new, more complex, descendant 
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Figure 8: Graphs Ti, T 2 and as described in Definition 14.31 



by replacing a cracker in the original descendant with two new crackers. We say that such a 
new descendant has been obtained from parthenogenesis. For simplicity of terminology, we again 
refer to the original descendant as the parent of the new descendant, and likewise we refer to the 
new descendant as the child of the original descendant. Also for simplicity of terminology, we 
refer to the the three breeding operations and the three parthenogenic operations collectively as 
the six breeding operations. 

Definition 4.4. A type 1 parthenogenic operation is a function V\ defined on the tuple (r l5 ei) 
where Ti = (Vi,Ei) is a bridge graph and e\ = (a, b) S E\ is a 1-cracker. This function maps 
such a tuple onto another tuple (r '.d, Vi,i>4) as follows 

where Td — (Vd,Ed) and {vi,v 4 } G Vd- The new set of vertices is Vd = V\ U {v\,V2,V3,Vi}. 
The new set of edges is E D = (Ei\ei) U {(a,vx), (vi,v 2 ), (vi,v 3 ), (v2,v 3 ), (v2,Va), (u3,«4), 
(vi,b)}. This process inserts an additional 1-cracker into Td- 




Figure 9: Graphs Ti and as described in Definition 14.41 



We refer to the subgraph T s = ( { «i , v 2 , v 3 , v 4 } , { («i , v 2 ) , («i , v 3 ) , (v 2 , v 3 ) , (« 2 , v 4 ) , (« 3 , «4 ) } ) as the 
parthenogenic diamond, and say that a type 1 parthenogenic operation inserts a parthenogenic 
diamond into a bridge. See Figure [S] for an illustration. 

Definition 4.5. A type 2 parthenogenic operation is a function V 2 defined on the tuple (ITi, ei, 62) 
where Y\ = (Vi, is a cubic graph containing a 2-cracker comprising two edges ei = (a, fe) and 
G2 = (c,d). This function maps such a tuple onto another tuple (Td,vi,V2) as follows 



^2(ri,ei,e 2 ) = (T D ,v 1 ,v 2 ), 
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Figure 10: A parthenogenic diamond 



where To — {Vd,Ed) and {vi,v 2 } £ Vd- The new set of vertices is Vd = V\ U {v\,v 2 }. The 
new set of edges is Ed = {E\\{e\, e 2 }) U {(a, V\), (b, (c, v 2 ), (d,v 2 ), (i>i) V2)}. This process 
inserts an additional 2-cracker into r^j- 




Figure 11: Graphs Ti and Td as described in Definition 14.51 



We refer to the subgraph T$ = ({vi, v 2 }, (v\, v 2 )) as the parthenogenic bridge, and say that a 
type 2 parthenogenic operation inserts a parthenogenic bridge into a 2-cracker. See Figure [11] 
for an illustration. 




Figure 12: A parthenogenic bridge 

Definition 4.6. A type 3 parthenogenic operation is a function V3 defined on the tuple (Ti,a) 
where Ti = (Vi, E\) is a bridge graph and a £ Vi is a vertex incident to a 1-cracker composing 
an edge e\ = (a, b) £ Ex and is adjacent to vertices c and d £ V\. This function maps such a 
tuple onto another tuple (To, a, Vx, v 2 ) as follows 

^3(^1, a) = (^D,a,n,n), 

where Td = (Vd,Ed) and {a, V\, v 2 } £ Vd- The new set of vertices is Vd = Vi U {vi,v 2 }. The 
new set of edges is E D = (Ex\{(a,c), (a,d)}) U {(o,ui), (a, u 2 ), (vx,v 2 ), («i,c), («2,d)}' 

We refer to the subgraph = ({a, i>i, U2}, ui), (a, V2), (wi, V2)}) as the parthenogenic trian- 
gle, and say that a type 3 parthenogenic operation inserts a parthenogenic triangle next to the 
1-cracker. See Figure [TBI for an illustration. 
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Figure 13: Graphs Ti and as described in Definition 14.61 




Figure 14: A parthenogenic triangle 



Lemma 4.7. A child graph resulting from any of the six breeding operations is connected. 

Proof. The nature of the six breeding operations is that the parent graphs are mostly unaltered, 
and are changed only in a neighbourhood of the introduced cubic cracker. Therefore, we can 
focus just on these areas. Since the parent graphs are (by definition) connected to begin with, 
we only need to be concerned with which edges, present in the parent graphs, are not present in 
the child graph. 

For the cases of type 2 breeding, and types 1 and 2 parthenogenesis, only a single edge from the 
parent graph (or from each of the parent graphs in the case of type 2 breeding) is missing in 
the child graph. By definition this edge cannot be a 1-cracker, and therefore, the parent graphs 
remain connected, and by construction it is clear that the adjoining cracker ensures the child 
graph is also connected. 

For type 1 breeding, only a single edge is removed from each parent graph. If neither edge is 
a 1-cracker, then the argument in the previous paragraph can be used to show the child graph 
is connected. However, it is possible that one or both removed edges could be 1-crackers. If so, 
the corresponding parent graphs become disconnected. However, if this is the case, the graphs 
are reconnected by the introduction of vertices V\ and V2 (see Definition 14. lj) . It is then clear by 
construction that the adjoining cracker ensures the child graph is also connected. 

For type 3 breeding, a vertex is removed from both parent graphs. From Definition 14. 3[ we 
know that none of the edges adjacent to these two vertices constitute 1-crackers. Therefore, the 
removal of these vertices cannot disconnect either graph. By construction it is then clear that 
the adjoining cracker ensures the child graph is also connected. 

Finally, for type 3 parthenogenesis, although two (adjacent) edges from the parent graph are 
missing in the child graph, it is clear from the latter's construction that this can not result in a 
disconnected descendant. □ 
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4.3 Inverse breeding and inverse parthenogenic operations 



For some tuples (r c ,e) where {e} G E D is a 1-cracker, the inverse operation B^ 1 (T D ,e) = 
(ri, T2, ei, e 2 ) is well defined. In such a case {e} is called an irreducible 1-cracker. If not 
{e} will be called a reducible 1-cracker. Similarly if the inverse operation B^^d, e 3 , e 4 ) — 
(Fx, r 2 , ei, e 2 ) is well defined, where {e 3 ,ei\ G Ed is a 2-cracker, the 2-cracker is called an 
irreducible 2-cracker and reducible 2-cracker otherwise. We will show later that the inverse 
operation B 3 1 (Td, ei, e 2 , 63) = (ri,r2, V\,V2, a, b, c, d, e, /) is always defined where {ei, e 2) 63} G 
Ed is a 3-cracker. Therefore every 3-cracker is irreducible. 

Definition 4.8. Whenever a cubic cracker is irreducible one of the equations (QJ, © and ([3]) 
defines the corresponding inverse breeding operation , B 2 1 (.) or B 3 The two cubic 

graphs (Ti,r2) from the tuple produced by these operations are parents of To- In particular, 

B^\T D ,e) = (T 1 ,T i ,e 1 ,e 2 ) (1) 

where Tp, e, Ti, T 2 , &i and e2 are defined in Definition 14.11 Similarly, 

B 2 ~ 1 (r D ,e 3 ,e 4 ) = (r 1 ,r 2 ,e 1 ,e 2 ) (2) 

where Tp, e^, e^, T%, r 2 , e% and e2 are defined in Definition 14.21 Also, 

#3 1 (Td, ei, e 2 , e 3 ) = (r 1 ,r 2 ,v 1 ,v 2 , a, b, c, d, e, /) (3) 

where Td, ei, e 2 , e 3 , Ti, T 2 , vi, v 2 , a, b, c, d, e and / are defined in Definition 14.31 

Similarly, inverse parthenogenic operations can be defined as follows. 

Definition 4.9. Equations (fj|, © and © define the corresponding inverse parthenogenic op- 
erations Vi 1 (.) , P 2 1 ( ) or V3 (.)■ The cubic graph (ri) from the tuple produced by these 
operations is a parent of Td- In particular, 

Vi 1 (T D ,v 1 ,v 4 ) = (T 1 ,e 1 ), (4) 

where Td, v\, v 4 , Ti and e± are defined in Definition 14.41 Similarly, 

V 2 1 {T D ,v 1 ,v 2 ) = (r 1 ,ei,e 2 ), (5) 
where To, vi, v 2 , Ti, ei and e 2 are defined in Definition 14.51 Also, 



V 3 1 (T D ,a,v 1 ,V2) = (Ti,a), (6) 
Where Td, a, Ti, T 2 , v% and v 2 arc defined in Definition 14.61 

Collectively, we refer to the three inverse breeding operations and the three inverse parthenogenic 
operations as the six inverse operations. 

It is important to note that the six breeding operations and the six inverse operations pre- 
sented here are not entirely new, and have been used in various forms in other cubic graph 
generation routines. For example, type 2 parthenogenesis induces an H-subgraph which is well- 
studied in literature (e.g. see Ore |Ore 67] ). Type 1 inverse parthenogenesis and type 1 in- 
verse breeding appear as Operation 02(^4) and Operations 1Z respectively in Ding and Kanno 
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|Ding and Kanno 06| . Types 1 and 2 parthenogenesis appear in Brinkmann [Brinkmann 96] . All 
of the six breeding operations except type 3 breeding appear in some sense as generating rules 
in Batagelj [Batagclj 86 , specifically generating rules PI, P2, P3, P4 an P8. However, all of 
the above works involve growing the complexity of a single graph by the evolution of subgraphs, 
rather than combining several cubic graphs together. In addition, the operations in the above 
works that are analogous to our parthenogenic operations are not confined to the same conditions 
as ours (that is, they must occur on 1-crackers and 2-crackers). Little consideration is given in 
the above works to procedures that are analogous to our inverse operations. The benefits and 
potency of the particular set of breeding and inverse operations that we have detailed in this 
section are demonstrated in the following section. 



5 Results 



The following three propositions relate to the different possible methods of creation of cubic 
crackers by the six breeding operations, and are used in the proof of the main theorem for this 
section, Theorem 15.51 

Proposition 5.1. Any descendant involving a 1-cracker can be obtained from either type 1 
breeding, type 1 parthenogenesis, or type 3 parthenogenesis. 

Proof. Consider a descendant Yd = (Vd,Ed) containing a 1-cracker comprising an edge 
(vi, v%). Since is cubic, v\ and V2 will both be adjacent to two more vertices, say {a, b} and 
{c,d} respectively. Since the 1-cracker comprises edge (vi,V2), we know that {a, b} and {c, d} 
are disjoint sets. Then, we can consider two cases. 

Case 1: The edges (a, b) and (c, d) are not present in To- In this case, the 1-cracker is irre- 
ducible. Suppose the bridge {v\,v%) is removed from Yd, separating the graph into two sub- 
graphs, Tsi — (Vsi,Esi) and Y$2 — (Vs2,Es2)- Without loss of generality, we assume that 
v\ G Vsi, and V2 £ Vs2- Then, we define a cubic graph Y\ = (Vi, E±), where V± = Vsi\vi and 
Ei = (Esi\{(a,vi), (b,vi)}) U (a.b). Similarly, we define a second cubic graph T 2 = (y2,E 2 ), 
where V2 — Vs2\v2 and E 2 = (i?s2\{(c, V2), (d, V2)}) U (c, d). Then, can be obtained from 
the type 1 breeding operation £>i(Ti,r2, (a, 6), (c, d)). 




Figure 15: Graphs Fx, T 2 and as described in Case 1 of Proposition l5.ll 



Case 2: At least one of the edges (a, 6) or (c, d) is present in Yd. In this case, the 1-cracker 
cannot be obtained from a type 1 breeding operation, as such an operation would remove these 
edges. A 1-cracker of this type is, therefore, reducible. If both edges are present, we can focus 
on either one. Without loss of generality, we will assume that edge (a, b) € Ed- Then, since Yd 
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is cubic, and vertices a and b are both adjacent to vertex v\ and to each other, they will also be 
adjacent to one more vertex each, say vertices e and / respectively. Note that is is possible that 
e = /, so we need to consider the cases separately. 

Case 2.1: e = /. Then, both edges (a, e) and (b, e) are in Ed, as seen in the right panel of 
Figure [TBI Since is cubic, vertex e is adjacent to a third vertex, say g. Clearly, edge (e,g) 
is a bridge, since edge (vi,i>2) is a bridge. Then, we define a cubic graph Ti — (Vi,Ei), where 
Ei = (E D \{(vi,v 2 ), (a,Vx), (Mi), ( a , e ), (M), (e,flO}) U (g,u 2 ) and V x = Vb\{vi, a, b, e}. 
We can then obtain Td from the type 1 parthenogenic operation ^(Ii, (<7, V2)). 




Figure 16: Graphs Ti and T^j as described in Case 2.1 of Proposition 15. II 



Case 2.2: e^/. Then, as illustrated in Figure [T71 we define a cubic graph T x = (Vi, Ex), where 
E 1 = (E D \{(a, vi), (b, Vl ), (a, 6), (a, e), (b, /)}) U {(e, (/, w x )} and Vi = 7 \{a, &}. Then, we 
can obtain from the type 3 parthenogenic operation T > a(Ti,Vi). □ 




Figure 17: Graphs Ti and as described in Case 2.2 of Proposition 15. II 



Proposition 5.2. Any descendant involving a 2-cracker can be obtained from either type 1 breed- 
ing, type 2 breeding, type 1 parthenogenesis, type 2 parthenogenesis, or type 3 parthenogenesis. 

Proof. Consider a descendant = (Vd, Ed) containing a 2-cracker comprising edges (v\, V2) 
and (w3,«4). We consider two cases. 

Case 1: Neither edge (uijWa) nor (v2,V4) are not present in Trj, as illustrated in Figure [18J In 
this case, the 2-cracker is irreducible. Suppose the edges {v\,V2) and (^3,^4) are removed from 
To, separating the graph into two subgraphs Tsi = (Vsi,Esi) and rs 2 = (Vs2,Es2)- Without 
loss of generality, we assume that vi,v 3 G Vsi and t>2,«4 6 Vs2- Then, we define a cubic graph 
Ti = (Vi,Ei), where V\ = Vsi and E\ — E$i U (vi,v 3 ). Similarly, we define a second cubic 
graph T 2 = (V2,i?2), where V2 = Vs2 and E 2 = Es2 U (i>2,W4). Then, Td is obtained from the 
type 2 breeding operation B 2 {T\, T 2 , (t>i, v 3 ), (v 2 , V4)). 
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Figure 18: Graphs Ti, and as described in Case 1 of Proposition 15.21 



Case 2: At least one of the edges (vi,v 3 ) or (112,^4) is present in Td- In this case, the 2-cracker 
cannot be obtained from a type 2 breeding operation, as such an operation would remove these 
edges. A 2-cracker of this type is therefore reducible. If both edges are present, we can focus on 
either one. Without loss of generality, we will assume that edge (vi,v 3 ) £ Ed- Then, since Yo 
is cubic, and vertex v± is adjacent to vertices vi and U3, it will be adjacent to one more vertex, 
say vertex a. Similarly, since vertex v 3 is adjacent to vertices v\ and m, it will be adjacent to 
one more vertex, say vertex b. Note that it is possible that a — b, so we need to consider the 
cases separately. 

Case 2.1: a = b. In this case, edge («i,a) G Ed and edge (v 3 ,a) E Ed, as illustrated in 
Figure [TSJ Since Td is cubic, vertex a is adjacent to a third vertex, say c. Clearly, edge 
(a, c) is a bridge. Then, we define a cubic graph Ti = (Vi,Ei), where V± = Vd\{vi,Vs} and 
Ei = E D \{(v 1 ,v 2 ),(v 3 ,v i ),(v 1 ,v 3 ),{a,v 1 ),(a,v 3 )} U {(a, v 2 ), (a, ^4)}- Then, we can obtain T D 
from the type 3 parthenogenic operation V 3 (Ti, a). Note that this case is essentially the same 
as Case 2.2 in Proposition [5TTJ 




Figure 19: Graphs fi and Td as described in Case 2.1 of Proposition 15.21 



Case 2.2: a ^ b. It is obvious that the non-adjacent edges (t>i,a) and (v 3 ,b) form a cutset. 
This implies that either both edges are 1-crackers, or together they form a 2-cracker. The former 
situation is covered in Proposition ?? and Td can be obtained from either type 1 breeding, type 
1 parthenogenesis, or type 3 parthenogenesis. For the latter situation, let us define a cubic graph 
ri = (Vl,Ei), where V\ = Vd\{wi,w 3 } and Ex = E D \{(vi,V 2 ), (^3,^4), (vi,v 3 ), (a,ui), (M3)} u 
{(a, V2), (b, V4)}, as illustrated in Figure [20l We can then obtain Td from the type 2 parthenogenic 
operation Vzfii, (a, w 2 ), (b, V4)). □ 

Proposition 5.3. Any descendant involving a 3-cracker can be obtained from type 3 breeding. 

Proof. Consider a descendant T d = {Vd, Ed) containing a 3-cracker comprising edges (i>i, W2), 
(v 3 ,Vi) and (v5,ve), as illustrated in Figure I2T1 Then, suppose the edges (fi,?^), {v 3 ,v<t) and 
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Figure 20: Graphs Ti and Yr> as described in Case 2.2 of Proposition 15.21 



(v$,vq) are removed from Tp, separating the graph into two subgraphs Y$i — (Vsi,Esi) and 
Ts2 = (Vs2, Eg2)- Without loss of generality, we assume that {v\, v 3 , v$} £ V51 and {v 2 , V4, v e } 6 
V51. Then, we introduce a new vertex V7, and dehne a cubic graph Ti — (Vi,E±), where 
V\ = Vsi U {^7} and E\ = Egi U {(vx, U7), (1)3, V7), (1)5, V7)}. Similarly, we introduce a 
new vertex wg, and define a second cubic graph T2 = (V^,^), where V2 — V52 U {v&} and 
E2 = Es2 U {(v2, Vg), («4, v%), (wg, v$)}. Then, we can obtain from the type 3 breeding 
operation B 3 (Tx, T 2 , v 7 , v 8 , vi, v 3 , w B) v 2 , V4, v e )- □ 




Figure 21: Graphs Ti, T2 and Td as described in Proposition 15.31 

Note that any 3-cracker can be obtained from a type 3 breeding operation, and therefore all 
3-cr ackers are irreducible. 

Definition 5.4. We refer to a set of genes that can, through a series of breeding and parthenogenic 
operations, be used to produce a descendant as ancestor genes for Tp- 

Propositions 15.11 - 15.31 and Definition 15.41 allow us to derive the main theorem of this section. 
Theorem 5.5. Consider any descendant cubic graph Tjj. Then, 

(1) Yd can be obtained from one or two parents by at least one of the six operations fii(-), 

b 2 (-), Bs(-), Vi{-), r 2 (-), 

(2) Every descendant Yd has a set of ancestor genes. 



Proof. Any descendant graph contains at least one cubic cracker. Any one of these cubic 
crackers can be selected, which will be either a 1-cracker, a 2-cracker, or a 3-cracker. Then, from 
Propositions 15.1] 15. 21 and I5.3[ we can obtain Yd, from either one or two parents, using one of the 
six operations. Therefore (1) is proved. 
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Next, consider how Yd is obtained. If Tjj is obtained through breeding, it has two parent graphs. 
If Td is obtained through parthenogenesis, it has one parent graph. These parent graphs may 
be either genes, or descendants. If any of them are descendants, then by part (1) they also 
have one or two parent graphs each. Inductively, we can continue to consider the parents of 
descendants, while recording which operations are used to produce them, thereby obtaining an 
ancestral family tree of Trj. Once the entire tree is determined, all of the top nodes are genes, 
and we can recall the sequence of operations that produces T]j from these genes. Therefore (2) 
is proved. □ 

See Figure [52] for an example of an ancestral family tree for a descendant with 14 vertices. 




Figure 22: An ancestral family tree for a 14- vertex descendant, with three ancestor genes 

Theorem 15.51 indicates that, for any descendant, we can obtain a set of ancestor genes by first 
applying an inverse operation to obtain one or two parents. Then we can apply an inverse 
operation on the parent (s) to obtain new parents (grandparents of the original descendant), and 
continue to apply inverse operations until a set of ancestor genes is obtained. 

However, a given descendant may contain only reducible cubic crackers, which do not permit 
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inverse breeding operations. In these cases, inverse parthenogenic operations can be performed, 
if one of the three parthenogenic objects are present within the descendant. The removal of such 
a parthenogenic object often changes a reducible cracker into an irreducible cracker which, in 
turn, permits an inverse breeding operation to be carried out. The following proposition ensures 
that, for any descendant, at least one of the cubic crackers permits an inverse operation, and 
furthermore that every reducible cracker in a given descendant can be changed into an irreducible 
cracker by a sequence of inverse parthenogenic operations. 

Proposition 5.6. If a graph To is a descendant, one of the following must be true. 

(1) To contains at least one irreducible cubic cracker, or 

(2) It is possible to perform a sequence of inverse parthenogenic operations, each time obtaining 
a new parent, until a parent is obtained that contains at least one irreducible cubic cracker. 



Proof. Since To is a descendant, it contains one or more cubic crackers. If any of them are 
irreducible, then (1) is true. If all the cubic crackers are reducible, then To contains at least 
one 1-cracker or 2-cracker (as 3-crackers are always irreducible). We can therefore select cither 
a reducible 1-cracker or a reducible 2-cracker in To- We will consider both cases separately. 

Case 1: We select a reducible 1-cracker C\ = {(a,b)}. Since To is cubic, we know that a 
must be adjacent to two more vertices, say c and d. Since C\ is reducible, then without loss of 
generality, we can assume that edge (c, d) is present in r^. Then, the cubicity of To also ensures 
that vertices c and d must each be adjacent to one more vertex, say e and / respectively. Note 
that it is possible that e = /. 

Case 1.1: If e = /, then the cubicity of To ensures that vertex e must be adjacent to one 
more vertex, say g (see Figure[23]). Then, removing edges (a, b) and (e,g) from To disconnects a 
parthenogenic diamond, which we can remove from by use of the inverse type 1 parthenogenic 
operation Vi 1 {Td, (a, b), (e, g)). In the parent graph Tp, there is a new 1-cracker comprising 
edge (b,g). Note that the new 1-cracker may also be reducible. 




Figure 23: Graphs and Tp as described in Case 1.1 



Case 1.2: lie ^ f, then removing edges (a, b), (c, e) and (d, f) from To isolates a parthenogenic 
triangle (see Figure which we can remove from by use of the inverse type 3 parthenogenic 
operation V% (Td> a, c, d). In the parent graph Tp, the original 1-cracker C\ = {(a, b)} remains. 
Note that the 1-cracker may still be reducible. 

Case 2: We select a reducible 2-cracker C2 — {(a, b), (c, d)}. Since C2 is reducible, then without 
loss of generality, we can assume that edge (a, c) is present in To- Then, the cubicity of To also 
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Figure 24: Graphs and Tp as described in Case 1.2 



ensures that vertices a and c must each be adjacent to one more vertex, say e and /, respectively. 
Note that it is possible that e = /. 

Case 2.1: If e = /, then as illustrated in Figure [25l the cubicity of Tp ensures that vertex e 
must be adjacent to one more vertex, say g. Then, removing edges (a, 6), (c,d) and (e,g) from 
Tp> disconnects a parthenogenic triangle, which we can remove from Td by use of the inverse 
type 1 parthenogenic operation 7- > 1 ~ 1 (Fd, a, c). In the parent graph Tp, a 1-cracker comprising 
edge (e, g) remains. Note that this 1-cracker may still be reducible. 




Figure 25: Graphs Yd and Tp as described in Case 2.1 

Case 2.2: If e ^ /, then as illustrated in Figure [2^1 the removal of edges (a, b), (c, d), (a, e) and 
(c, /) from disconnects a parthenogenic bridge, which we can remove from by use of the 
inverse type 2 parthenogenic operation (T]j,a, c). In the parent graph Tp, there is a new 
2-cracker comprising edges (e,b) and (f,d). Note that this 2-cracker may be reducible. 




Figure 26: Graphs Tp> and Tp as described in Case 2.2 

Since in all cases considered, the parent graph contains a cubic cracker (either a new one intro- 
duced by an inverse parthenogenic operation, or one that remains from the original descendant), 
the parent is itself always a descendant. Then, either (1) is true for this parent, or if not, we 
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can repeat the above procedure until we obtain a parent for which (1) is true. Since each inverse 
parthenogenic process outputs a parent with fewer vertices than its child, we are guaranteed to 
eventually converge to such a case. □ 



6 Conclusions 



The theory presented above allows us to separate the set S of connected cubic graphs into two 
distinct and encompassing categories - the comparatively smaller set of genes, that form the 
basic building blocks of S, and the much larger set of descendants, that inherit a lot of their 
structure from the genes. Theorem 15.51 and Proposition [5]6] give both a proof of existence, and a 
guarantee of obtaining a set of ancestor genes for any given descendant. An algorithm to identify 
a set of ancestor genes, given a descendant, would be a simple task of identifying all the crackers, 
surveying each until one is found that permits an inverse operation, and recursively repeating 
the process in each parent obtained until only genes remain. Such an algorithm would terminate 
in polynomial time. 

In addition to the aesthetic beauty of rendering the ancestry of cubic graphs as finite sets of 
smaller cubic graphs, there are some important algorithmic benefits as well. Since descendants 
inherit much of their structure from genes, it is possible that a search for graph theoretic prop- 
erties within a descendant could be more efficiently recast as a search of the (typically much 
smaller) genes instead. In this context, Theorem 15.51 and Proposition 15.61 give rise to a generic 
decomposition algorithm that could be applied on most cubic graphs. Since the ancestor genes 
also lie within the set of connected cubic graphs, any existing algorithms designed for cubic 
graphs will work for the genes. 

For example, experimental evidence indicates that many non-bridge non-Hamiltonian cubic 
graphs are descendants that contain at least one mutant ancestor gene. In these cases, it is 
clear that the descendant has inherited the non-Hamiltonicity property from its ancestor mutant 
gene (or genes). Then, an obvious heuristic for determining Hamiltonicity in a descendant is to 
identify a set of ancestor genes and determine their Hamiltonicity instead, using whatever state 
of the art algorithms are available (e.g. see Eppstein Eppstcin 03 ). Since determining Hamil- 
tonicity is an NP-complete problem, and therefore the best known algorithms have exponential 
solving time, such a decomposition represents a large saving in solving time. It is potentially 
possible that the Hamiltonian cycle problem, already known to be NP-complete even when con- 
sidering only cubic graphs [Garey and Johnson 79| , could be further refined to requiring the 
consideration of only genes. If so, the extreme rarity of mutants indicates such an avenue could 
potentially be fruitful. The inheritance of non-Hamiltonicity, and other such graph theoretic 
properties, is a subject for future research. 

Given that properties can be inherited from a set of ancestor genes, a natural question to ask 
is how many different sets of ancestor genes a descendant graph might have, and how their 
various (potentially conflicting) properties may influence the single descendant. The following 
conjecture, if correct, removes any such confusion. 

Conjecture 6.1. Any descendant Td has a unique set of ancestor genes. 

To support Conjecture 16.11 we conducted experiments on individual graphs in which we con- 
sidered all possible ways to decompose the graph into a set of ancestor genes, and verified that 
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each of these approaches gave the same set of ancestor genes. This experiment was performed 
on all cubic graphs containing up to 18 vertices, and all triangle-free 20-vertex cubic graphs, 
constituting 143,528 graphs. For these sizes, there can be up to 7! = 5040 possible ways to 
decompose a graph into ancestor genes. No counterexample to Conjecture 16. II was found among 
the tested graphs. 

If Conjecture 16.11 is true it provides, along with Theorem 15.51 a guarantee of existence and 
uniqueness of ancestor genes for every cubic graph. For such a graph, the above implies that it 
is either a gene, or there is a unique set of ancestor genes which can be identified in polynomial 
time, and that these ancestor genes provide the majority of structure in the descendant. If 
false, Conjecture 16.11 may still be true for all but very special classes of descendants. Note that 
although the conjecture postulates the existence of a unique set of ancestor genes for any given 
descendant, the order of breeding operations used to obtain the descendant is clearly not unique. 
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